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Thermodynamics of an attractive 2D Fermi gas 
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Thermodynamic properties of matter are conveniently expressed as functional relations between 
variables known as equations of state. Here we experimentally determine the compressibility, density 
and pressure equations of state for an attractive 2D Fermi gas in the normal phase as a function 
of temperature and interaction strength. In 2D, interacting gases exhibit qualitatively different 
features to those found in 3D. This is evident in the normalized density equation of state, which 
peaks at intermediate densities corresponding to the crossover from classical to quantum behaviour. 
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Two-dimensional (2D) quantum matter can display be¬ 
haviors not encountered in three-dimensions (3D) ^[2]. 
Fermions confined to 2D planes play a key role in uncon¬ 
ventional superconductors [3], graphene [1] and certain 
topological insulators [S] yet understanding the proper¬ 
ties of these complex materials can present significant 
theoretical challenges. Simpler systems, such as ultra¬ 
cold atomic gases with tunable interactions [BHH] , can 
serve as valuable testbeds for building up and validat¬ 
ing models of interacting fermions in 2D. Experiments 
on 2D Fermi gases to date have addressed their produc¬ 
tion isHn], pairing nmn], pseudogap m and polaron 
physics [151 [H] j pair condensation the Berezinskii- 
Kosterlitz-Thouless transition [18] and the pressure [19] 
as a function of interaction strength at low temperatures. 
Theoretically, several groups have investigated superflu¬ 
idity [3DH33] and the thermodynamic properties of 2D 
Fermi gases [3H1 - I37] , however, a full experimental charac¬ 
terization remains to be established. 

In cold atom experiments, a 2D gas can be produced 
by subjecting a cloud to tight transverse confinement. 
In a harmonic potential, with characteristic frequencies 
u}x,ujy and a gas can become kinematically 2D when 
the transverse (z) confinement energy exceeds all other 
energy scales including the thermal energy, chemical po¬ 
tential and interaction energies. The first two criteria 
require that huiz ^ kBT,EF where h is Planck’s con¬ 
stant, ks is Boltzmann’s constant, T is the temperature 
and Ep is the Fermi energy. The criterion for the in¬ 
teractions requires that neither elastic or inelastic (e.g. 
pair-forming) collisions result in the population of trans¬ 
verse excited states [SlEsHSg. At nonzero densities it 
was found that moderate interactions can lead to trans¬ 
verse excitations, even when kpT and Ep lie below the 
transverse oscillator energy [33j . This in turn can im¬ 
pact other parameters such as the pairing gap |34j and 
superfluid transition temperature [24j . 

Here we measure the thermodynamic properties of at¬ 
tractive Fermi gases in the normal phase in a parameter 
regime where 2D kinematics has been clearly established 
[33] . We apply a protocol based on the approach used for 


a 3D Fermi gas at unitarity |35j and employed on a 2D 
Bose gas [33]. By establishing a model independent re¬ 
lationship between the compressibility and pressure, we 
extract a range of thermodynamic properties including 
the temperature, chemical potential and density equa¬ 
tions of state. 

In a two-component Fermi gas thermodynamic vari¬ 
ables can be expressed as functions, /, connecting the 
underlying energies within the system m- At fixed tem¬ 
perature and interactions, these can be related to the 
density via the Gibbs-Duhem equation. In 2D the pres¬ 
sure P, density n and isothermal compressibility k are 
given by 


P = = J T)dn' 
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where fi{l3pL, /SEb) depend on the chemical potential /i, 
temperature and interaction strength, /3 = \/{kBT), 
A = I(mkBT) is the thermal de Broglie wave¬ 

length, m is the atomic mass and Eb is the two-body 
binding energy which in quasi-2D is governed by coz and 
the 3D scattering length a. The 2D scattering length a 2 D 
is related to the binding energy by Eb = h^/{ma 2 p,) in 
the kinematically 2D regime [6l HI [28] . Knowledge of the 
functions f; represents a complete determination of the 
thermodynamics and can establish valuable benchmarks 
for comparison with many-body theories |38j . 

In the experiments that follow we study 2D atom 
clouds at thermal equilibrium in a cylindrically sym¬ 
metric harmonic trap Vr(x,y) = 'muy^{x^ y^)l‘2 with 
ujr{= lOx = oJy) ujz. Due to the slowly varying ra¬ 
dial potential we can make use of the local density ap¬ 
proximation (LDA) which asserts that local thermody¬ 
namic properties will be equivalent to those of a homo¬ 
geneous gas at the same temperature and chemical po¬ 
tential [33]. Under the LDA the atomic density can be 







2 



Figure 1. (a) Average of 10 absorption images of a 2D Fermi 
gas, with l3Eb = 0.26, prepared under the same experimental 
conditions {B — 880 G, N = 16,000, T « 25 nK) (b) Az- 
imuthally averaged density, n(14), obtained from image (a) 
as a function of the local potential, (c) Dimensionless com¬ 
pressibility K vs. dimensionless pressure p of a 2D Fermi gas 
with PEb = 0.26(2). Faint green circles show data extracted 
from single images and dark green triangles are data binned 
according to p. Grey solid line shows the equation of state 
for a noninteracting Fermi gas and dashed green line is the 
predicted k based on the virial expansion to third order. 


written as n[/x(a;,y),T], where p,{x,y) = /xq — Vr{x,y) 
and pq is the chemical potential at the trap center. In 
any single realisation of the experiment, the parameters 
P and Eb will be fixed across the cloud such that Eqs. (1- 
3) yield fn{Py,pEb) = fp{py,PEb) and f^{Py,PEb) = 
fp{Pfi, PEb) by differentiation with respect to Py. 

We prepare single 2D clouds of neutral ®Li atoms in 
a hybrid optical/magnetic trap at temperatures in the 
range of 20-60 nK |33l |40H42]. A blue-detuned TEMqi 
mode laser beam provides tight confinement along z with 
a;z/(27r) = 5.15kHz. Radial confinement arises from 
residual magnetic field curvature present when the Fes- 
hbach magnetic field is applied, leading to a highly har¬ 
monic and radially symmetric potential with lOrjpl'K') = 
26 Hz with an anisotropy of less than 0.6%. We image 
along z to directly obtain the density n{x, y) of clouds 
prepared in the kinematically 2D regime. Fig. 1(a) shows 
the average of 10 images taken under identical experi¬ 
mental conditions. Due to the symmetry of Vr{x, y), and 
a precise calibration of the absorption imaging |42j . we 
can azimuthally average these images to obtain n{Vr) as 
shown in Fig. 1(b). 

Both a and Vr{x,y) are precisely known for a given 
experimental sequence providing accurate knowledge of 
Eb and, via the LDA, the change in chemical potential 
dy{= —dVr) across the cloud. Other parameters how¬ 
ever, including the absolute temperature and chemical 





Figure 2. Normalized density equation of state for a 2D at¬ 
tractive Fermi gas for PEb = 0.47(3) (purple circles), PEb = 
0.26(2) (green triangles), PEb = 0.06(1) (blue daimonds) and 
PEb = 0.005(1) (red squares). Dashed lines show the calcu¬ 
lated equation of state using the third order virial expansion 
[49j . Solid blue line and grey stars are the calculated equation 
of state for PEb = 0.5 based on Luttinger-Ward (LW) [26] and 
quantum Monte-Carlo calculations (QMC) [27], respectively. 


potential, are unknown. Furthermore, absorption imag¬ 
ing is susceptible to systematics that make precise cali¬ 
bration of the atom density challenging [iHl - HH] . To pro¬ 
ceed we first obtain an estimate of the absolute temper¬ 
ature and chemical potential by fitting the low density 
wings of the cloud with the virial expansion in 2D to 
third order [15]. As only a small fraction of the cloud 
can be used, this can lead to relatively large uncertainties 
in the fits [42]. Fortunately, the relationships connect¬ 
ing fp, fn and /k allow this to be improved. In analogy 
with refs [351 ISSj, we use the n(lA) data to construct a 
model independent equation of state for the dimension¬ 
less compressibility k = k/kq versus pressure p = P/Pq 
[ 35] where kq = \/{nEF) and Pq = nEp/2 are the lo¬ 
cal compressibility and pressure of an ideal Fermi gas at 
T = 0, respectively, Ep = {fpTT/m)n = kpTp is the lo¬ 
cal Fermi energy and Tp is the Fermi temperature. In 
Fig. 1(c) we plot k against p for PEb = 0.26(2). At high 
temperatures the compressibility lies close to that for the 
ideal gas yet it deviates above this at lower temperatures 
(lower p). The virial expansion provides a reliable predic¬ 
tion of k for p> 6. Eqs. (l)-(3), along with dp = —dV^, 
then allow one to find the relative temperature T = T/Tp 
and /3/x at any position in the cloud using the integrals 


T = Ti exp 



(4) 


Pp = {Pp), 



(5) 


where the initial points can be chosen to lie in the range 
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where the virial expansion is accurate. Implementation 
of Eq. (4) turns out to be relatively insensitive to the pre¬ 
cise starting conditions and provides highly robust ther¬ 
mometry directly from the k vs. p equation of state. As 
a validation of the imaging calibration, the absolute tem¬ 
perature found from Eq. (4) should be consistent across 
the entire cloud and should match the value of that 
gave the best fit in the cloud wings using the virial ex¬ 
pansion [32] ■ 

The integral for the chemical potential, Eq. (5), should 
also yield values of /3/r that are consistent with temper¬ 
ature found using Eq. (4). Additionally, p should vary 
according to the LDA as a function of position across the 
cloud. Meeting all of these conditions requires accurate 
calibration of the absorption imaging as full consistency 
is only obtained with the correct parameters [42] . Agree¬ 
ment with the virial expansion at high temperatures pro¬ 
vides a complementary validation of the thermodynamic 
parameters that assures the data is free from systematics 
within the error bounds of the virial fit. 

With T/Tp and jdp known across the cloud, and the 
above criteria satisfied, we can construct the homoge¬ 
neous density equation of state. In Fig. 2 we plot n{p, T) 
for four values of pEh, normalized to the density of an 
ideal Fermi gas at the same chemical potential and tem¬ 
perature no{p,T) = (2/A^) ln(I-I-e^^). Reducing di¬ 
mensionality leads to a dramatic difference in this equa¬ 
tion of state compared to the 3D attractive Fermi gas 
[551150] . The normalized density n/no displays a non¬ 
monotonic behaviour as a function of /3/r. This is most 
evident for the gas with strongest interactions f3Ei, = 0.47 
where n/no peaks just below /3/r = 1. Also shown are 
the calculated equations of state using the virial expan¬ 
sion (dashed lines) [35], which is valid for /3p < —1.5, 
along with recent Luttinger-Ward (LW, solid line) [55] 
and quantum Monte-Carlo (QMC, grey stars) [57] cal¬ 
culations for a cloud with = 0.5. Our data for 
j3Eb = 0.47 shows good qualitative agreement with these 
calculations lying in between the LW and QMC curves. 

The peak in n/no originates from the presence of a 
two-body bound state which exists for arbitrarily weak 
attraction in 2D. Interactions will generally be most sig¬ 
nificant when the kinetic energy of colliding atoms, 
is approximately equal to Et- In 2D, when T ^ Tp, this 
occurs when the interaction parameter ln(kpa 2 D) 0 
[5], where kp = '/2jrn is the Fermi wavevector. At T = 0, 
the existence of a bound state leads to the possibility of 
tuning from the fermionic (BCS) regime where Eh <C Ep, 
to the bosonic regime with Eh ^ Ep, simply by varying 
the density. At finite temperatures, the thermal energy 
kpT sets a lower bound on the collision energy and the 
Bose limit is only possible when /iEh > 1. For the data 
in Fig. 2, PEh is always less than unity and Ej. always 
exceeds Eh- In the low density (classical) region of the 
clouds, where /3/r < 0, Ek will be set by kpT, and inter¬ 
actions become stronger as the density begins increasing 


and the relative temperature T decreases. However, once 
Pp > 1, Ep becomes the dominant energy (quantum 
regime) and E^ increases above kpT, and hence further 
above Eh- This leads to effectively weaker interactions at 
high density and n/no begins decreasing for pp > I. 



Figure 3. (a) Normalized free energy F and (b) internal en¬ 
ergy [/ as a function of T/Tp for PEh = 0.47(3) (purple cir¬ 
cles), PEh = 0.26(2) (green triangles), PEh = 0.06(1) (blue 
daimonds) and PEb — 0.005(1) (red squares). Inset in (a) 
shows a zoomed in view of F at low temperature. Dashed 
lines show the calculated energies using the virial expansion 
for the same values of PEb and grey lines show the ideal gas 
result, (c) Normalized pressure p as a function of the interac¬ 
tion parameter In {kpa^o) and temperature T /Tp. Solid lines 
show how the interaction parameter varies as a function of the 
relative temperature within single clouds and grey circles on 
the contours indicate the approximate location of the peak in 
n{p,T)/no{fJ.,T). 

Not all thermodynamic properties of an interacting 2D 
Fermi gas can be obtained from measurements on a sin¬ 
gle cloud. While the free energy, E = U — TS, where U is 
the internal energy and S is the entropy, is readily found 
via the grand potential fl = —PA = F — pN, where A 
is the area and N is the particle number, U and S can¬ 
not be found individually without further information. 
The normalized free energy is given by F/{NEp) = F = 
p/Ep —p/2, and is plotted in Fig. 3(a). However, unique 
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determination of U and S requires differentiation across 
clouds with different 020 ■ Specifically, by considering F 
as a function of temperature and a 2 D we can evaluate 
the contact density C using dF/d{lna 2 D) = ‘^C/kp at 
fixed T/Tp. With this the internal energy, and hence 
the entropy, can be found via the Tan pressure relation 
[KTVtHH] . p = 211/(NEp) + 2C/kp [32]. As we only have 
measurements at four values of a 2 D the measured con¬ 
tact contains significant uncertainty compared to other 
variables. However, as C remains relatively small in our 
experiments, the internal energy shows the expected be¬ 
haviour, Fig. 3(b). 

Finally, we can plot the dimensionless pressure p versus 
interactions and temperature in Fig. 3(c). Constructing 
a full 3D surface plot of p as a function of In {kpa 2 D) and 
T/Tp would represent a complete characterisation of the 
thermodynamics of the 2D Fermi gas. The contours in 
the p = 0 plane show how the relative temperature and 
interaction strength evolve in a single cloud. Grey circles 
on the contours indicate the approximate peak location 
in the density equation of state. 

In summary we have measured the thermodynamic 
properties of a 2D Fermi gas with attractive interactions 
in the normal phase. The existence of a bound state leads 
to qualitatively different behaviour to 3D gases as appar¬ 
ent in the density equation of state. For the gas with 
strongest interactions (/3Ab = 0.47) the superfluid tran¬ 
sition temperature is expected to lie at approximately 
0.05 T/Tf [2^ which is around a factor of two colder 
than we currently achieve. Future studies investigat¬ 
ing thermodynamic signatures of the superfluid transi¬ 
tion may provide insight into the nature of the transition 
in a quasi-2D Fermi gas. At stronger interactions one 
could investigate the effects of transverse excitations on 
the thermodynamic properties. 

We thank P. Hannaford, J. Drut, M. Parish and J. 
Levinsen for helpful discussions and comments on the 
manuscript and the authors of refs |2ni|27| for sharing 
their data. C. J. V. and P. D. acknowledge financial 
support from the Australian Research Council programs 
FT120100034, DP130101807 and DE140100647. 

Note added: A related work has recently been posted 
which examines the thermodynamic equation of state 
across the 2D BEC-BCS crossover [53] . 
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Preparing 2D Fermi gases 

To produce a 2D Fermi gas we begin with a 3D gas 
containing N/2 = = 2 x 10® Li® atoms confined in a 

red-detuned (1075 nm) single beam optical dipole trap. 
The cloud is evaporatively cooled to a temperature of 
T ~ 0.1 Tp at unitarity (832.2 G). We then ramp the 
magnetic field to the BEG side of the Feshbach resonance 
(780 G) and transfer the atoms to a TEMqi mode optical 
trap produced by a 532 nm frequency-doubled fibre laser 
[T]. This provides tight confinement in the transverse di¬ 
rection, while radial confinement is provided by residual 
curvature in the magnetic field generated by the Feshbach 
coils. Initially, the TEMqi mode laser is ramped up to 
approximately 10% of its full power over 350 ms and the 
3D optical trap is simultaneously switched off. We then 
perform evaporative cooling in the weak quasi-2D trap 
by applying a transverse magnetic field gradient over a 
period of 2 seconds. Including this extra step makes it 
possible to achieve lower temperatures and provides pre¬ 
cise control over the atom number in the final 2D trap. 
Lastly, the magnetic field gradient is removed and the 
TEMqi mode is ramped to full power to achieve a tight 
trapping frequency of 5.15 ± 0.24 kHz and a radial trap¬ 
ping frequency of 26.4 ±0.1 Hz at a magnetic field of 
972 G. After an equilibration time of 250 ms an in-situ 
absorption image is acquired. 


Calibration of the absorption imaging 


Absorption imaging of the trapped clouds is achieved 
by illuminating the atoms with a short pulse of resonant 
laser light propagating along the z-direction. The atoms 
leave a shadow in the laser light which is imaged onto a 
CCD camera. The atomic density in the x-y plane can 
be inferred from the resulting absorption image via the 
Beer-Lambert formula including saturation [5] 


n{x,y) 


1 

(7 



' If{x,yy 

.Io{x,y)_ 


+ 


h{x,y) - If{x,y) 


(SI) 


where If{x,y) is the measured intensity with the atoms 
present, Io{x, y) is the measured intensity in the absence 
of atoms, a is the absorption cross-section and Isat is the 
saturation intensity of the transition. Images used for the 
thermodynamic analysis are taken with a pulse length of 
10 ys and intensity I ^ 0.25 I^at- However, with these 
imaging parameters, the measured density generally un¬ 
derestimates the true density due to the Doppler shift 


of atoms accelerated by the imaging light and possible 
multiple scattering at finite gas density [3j. 

To correct for these effects we acquire a series of ad¬ 
ditional images of the 2D clouds using a 1 /xs imaging 
pulse with an intensity I ~ lOIsat- In this high intensity 
regime the second term in equation (SI I dominates and 
any Doppler shifts during this short pulse are negligible 
compared to the power broadened width of the transition. 
The measured optical density (OD) under these condi¬ 
tions provides a robust measure of the absolute density 
but contains significantly more statistical noise. By 
averaging many images taken in both low and high inten¬ 
sity regimes of clouds prepared in an identical way we can 
determine an overall correction to the optical density by 
constructing a correction factor Cqd = ODhigh/ODiow 
In general, this correction will be a function of the mea¬ 
sured optical density, however, as we only ever deal with 
very low optical densities (typically less than 0.1), Cqd 
turns out to be effectively constant over the range we 
access. In Eig. SI we plot the ratio ODhigh/ODiow ver¬ 
sus the measured ODiow in the region of the image where 
atoms are present. As can be seen the correction factor is 
approximately constant in the range of optical densities 
measured and becomes noisy at low OD (below ~ 0.02). 



ODiow 


Figure SI. Measured optical density correction factor Cod 
determined at a magnetic field of 972 G. Error bars represent 
the standard error determined for each bin of ODiow 

Averaging the measured correction factor in the range 
0.02 < ODiow < 0.06 and applying this to the images 
provides a corrected density distribution where the un¬ 
certainty in the corrected density is typically of order 
~ 10%. The actual correction factors determined in this 
way for the four magnetic fields investigated in the pa¬ 
per are shown in the central column of Table In the 
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next section, we show how fitting the density using the 
virial expansion provides a way to improve the accuracy 
of these coefficients. 


Temperature estimation using the 2D virial 
expansion 


The thermodynamic properties of an interacting 2D 
Fermi gas at relatively high temperatures T/Tp > 1, 
where Tp ]s the Fermi temperature, are well approxi¬ 
mated using the virial expansion [ 3 . When the fugacity 
j = is a small parameter, one can expand the grand 
potential in powers of the fugacity, where the expansion 
coefficient for the n-th order term depends on the solu¬ 
tion of the n-body problem for the interacting particles. 
For a 2D gas in a harmonic trap the density will be given 
to third order by 

n(r) = ^ (In [1 -H 5(r)] -h 2A&25(r)^ -f 3A635(r)^ H-) 

(S2) 


where A = 2tt h,/[inkpT) is the thermal de Broglie 

wavelength, m is the mass of the atom, j(r) = 
is the local fugacity, ^0 is the chemical potential at the 
centre of the cloud, V{r) = muj‘^{x‘^ + y‘^)j2 is the trap¬ 
ping potential, /3 = l/^kpT), tUr is the radial harmonic 
confinement frequency and A 62 and Ab^ are the 2 nd and 
3rd virial coefficients, respectively. The coefficients Ab^ 
are functions of pEb where Eb is the two-body binding 
energy that depends upon the transverse confinement fre¬ 
quency (jJz and the 3D scattering length a 0 H] . Eb is 
generally known precisely as both the trapping frequency 
and 3D scattering length are readily determined. 

The second order virial coefficient A 62 can be conve¬ 
niently calculated according to the Beth-Uhlenbeck for¬ 
malism. By introducing a dimensionless parameter 


a: = In 


A 

0-20 


= - In {2TTl3Eb), 


(S3) 


we obtain 


Ab2 = e^^'’- / dt 


exp 


27r 


TT^ + 4 (i — x) 


2 • 


(S4) 


The determination of the third order virial coefficient 
Ab^ is much more involved. It was calculated in Ref. [S] 
with the help of an isotropic harmonic confinement and 
in Ref. [5] by generalizing the diagrammatic approach in¬ 
troduced by Leyronas [S]. In Fig. S2, we show the numer¬ 
ical result. Empirically, in the range 0.003 < j3Eb < 1, 
we may use the fit. 


7 

A 63 = - ^ a„x", (S5) 

n=0 





Figure S2. Third order virial coefficient of a two-dimensional 
Fermi gas. 


Magnetic field (G) 

Cod (high/low) 

Cod (virial) 

972 

1.27 ±0.08 

1.21 ±0.03 

920 

1.24 ±0.17 

1.20 ±0.04 

880 

1.20 ±0.12 

1.21 ±0.03 

865 

1.24 ±0.13 

1.22 ±0.02 


Table I. Correction factors applied to the measured optical 
density of the 2D clouds at each magnetic field used in the 
experiment. The first estimate of the correction factor is given 
by Cod (high/low) which is determined by the ratio of the 
optical densities in both high intensity and low intensity im¬ 
ages. The correction factor listed is the average value of Cod 
in the range and the uncertainty is given by the standard er¬ 
ror. The values listed in column Cod (virial) give the scaling 
factor which minimises the difference between the measured 
density and density determined by fits to the virial expansion. 
Error estimates for these measurements is determined by the 
95% confidence interval of the fit. 


where Og = 0.45938, Oi = 0.40400, 02 = 0.31103, 
03 = 0.16998, 04 = 0.17801, 05 = 0.23461, og = 0.13623 
and 07 = 0.02685. With the ability to calculate these 
coeffecients we can fit equation S2 to the (azimuthally 
averaged) density profile with T and /ig as free parame¬ 
ters. 


The fitting procedure is sensitive to the absolute tem¬ 
perature T via the 1/A^ term and also through /?£!{, de¬ 
pendence of the virial coefficients. Additionally, any er¬ 
rors in the correction factor Cqd can lead to systematic 
errors in the fit. To account for each of these effects, we 
perform the fitting procedure iteratively using a bisec¬ 
tion algorithm for different values of (3Eb until the value 
of PEb used to determine the virial coefficients converges 
with the fitted temperature T. We also perform this fit¬ 
ting procedure with different values of the optical density 
correction factor Cqd to determine the combination of 
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PEb, T and Cqd which gives the best fit to the mea¬ 
sured density. This allows us to further refine the value 
of Cod as the quality of the fit is quite sensitive to the 
overall scaling factor. The right hand column of table 
|T] shows the optimised values of Co d that give the best 
fits at each magnetic field. The uncertainty in Cod is 
reduced to the few percent level by this procedure. In 
Fig. |S3[ we show an example of the resulting fit of the 
virial expansion to the cloud wings for PEj, = 0.005(1). 



Figure S3. Fit of the in-trap density with the third-order virial 
expansion. Blue points are the average of 20 experimental im¬ 
ages. Red solid line is the virial expansion fit. After conver¬ 
gence of the bisection algorithm the value of PEi, = 0.005(1) 
is found at an optimised Cod of 1 . 21 ( 2 ). 

At this point we have a first estimate of the abso¬ 
lute temperature and chemical potential of the 2D cloud 
which provides a basis upon which to build up the full 
thermodynamic description. 


Full thermodynamic analysis and validation of virial 

fits 


Equations (4) and (5) in the main paper provide a 
general approach to evaluate the relative temperature 
T = T/Tp and chemical potential /3/r via integration 
of the k vs. p equation of state at any temperature or 
interaction strength. Having performed the virial fit as 
described above, we can use the virial expansion to cal¬ 
culate k and p at the PEf, determined from the fit using 


K = 


2 /? 

{nXf 


-h 4A625^ + 9A535^ -I- 

1 + 5 


(S 6 ) 



oo 

J dt In (l -h je”*) -I- A 625 ^ -h Afesj^ -|- 


(S7) 

While these are valid only in the low density wings of the 
cloud where the relative temperature is high, they pro¬ 
vide useful initial conditions for evaluating the definite 


integrals for T and /3/x, equations (4) and (5). Rearrang¬ 
ing the expression for the density, equation (S2), we find 
that T is given by 


f = [in (1 + 5 ) + 2A625^ + 3A&35^ -I-] 


-1 


(S8) 


Under the assumptions of thermal equilibrium and the 
Local Density Approximation (LDA) we require that the 
absolute temperature across the cloud be uniform and 
that the chemical potential evaluated from equation (5) 
should satisfy /r(r) = /xq — U(r). To test compliance with 
these criteria, we first use the temperature integral, equa¬ 
tion (4), to find the relative temperature T for each value 
of p in the k-p equation of state. From this we can find 
the absolute temperature T using the local Fermi tem¬ 
perature Tp{r) = Ep{v)/kB, with Ep{y) = {nhP/m)n{r) 
being set by the density. Figure [S^a) shows the absolute 
temperature T determined at various positions through 
the cloud by varying the endpoint of the integration equa¬ 
tion (4). The data are approximately uniform across 
the entire cloud with some fluctuations appearing due to 
noise in the density profile n{r). The absolute tempera¬ 
ture is consistent with the virial result, grey line, where 
the uncertainty in the virial fit is indicated by the light 


grey band. Averaging the values in Fig. S4’a) over the 


different densities provides a robust estimate of the ab¬ 
solute temperature with a smaller uncertainty than that 
resulting from the virial fit. These provide the final val¬ 
ues of PEb which are given for each magnetic field in 
table mi 

A similar check can be performed to ensure the consis¬ 
tency of the chemical potential. It follows from the LDA 
that one can determine /xq from the value of /x(r) found at 
any position through the cloud using equation (5) from 
the main text. Again the virial expansion provides the 
initial values {PfJ.)i and Ti from which Pfi at any point 
in the cloud can be determined. As /3/x is calculated for 
each p we can infer the chemical potential in the trap 
centre using the temperature determined above and our 
knowledge of the trapping potential. Plotting the /Xq de¬ 
termined in this way allows us to validate both the LDA 
and the thermometry. Figure [S4j)b) shows a plot of /xg 
found by terminating the integral in equation (5) at dif¬ 
ferent positions through the cloud. The overall flatness 
of the curve confirms the consistency of the parameters 
found and used in the analysis. The fact that our results 
are fully consistent with the virial expansion in the high 
temperature limit provides a rigorous confirmation of the 
validity of the approach. 


Additional thermodynamic parameters 

Having determined pfj, and T/Tp, we can take the 
product of these to find the dimensionless chemical po- 
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Figure S4. Result of the temperature and chemical potential validation checks for a cloud with fiEb = 0.005(1). Panel (a) 
shows the resulting bulk temperature T calculated from different points in the cloud, whereas panel (b) is the calculated central 
chemical potential fio- Solid grey lines represents the value of T and /ro found via the fit to the wings of the in-trap density 
using the virial expansion and the shaded grey bands are the 95% confidence intervals of the fit. 


Magnetic field (G) 

PEb (virial) 

PEb (full analysis) 

972 

0.005(1) 

0.005(1) 

920 

0.06(1) 

0.06(1) 

880 

0.29(5) 

0.26(2) 

865 

0.49(5) 

0.47(3) 


Table II. Values of l3Eb determined at each magnetic field. 
PEh (virial) is found by via the fitting to the in-trap density 
distribution with the third order virial expansion after find¬ 
ing the optimal Cod- The quoted uncertainties represent the 
95% confidence interval of the fit. We also show the value of 
PEb as determined via the full thermodynamic analysis using 
the dimensionless compressibility and pressure. In this case 
the error includes both statistical uncertainty and systematic 
uncertainty arising due to measured field strengths and trap¬ 
ping frequency. 


tential ^/Ep which, in the case of a homogenous gas, is 
equivalent to the dimensionless Gibbs free energy. This 
is plotted in Fig. S5 for the four different interaction 
strengths. As can be seen, the chemical potential for 
interacting gases remains monotonic in the range of the 
normal phase that we measure but lies well below the 
ideal gas result even for modest interactions. 

Unlike the 3D Fermi gas at unitarity m, where the 
pressure follows a simple relation P — 2/3 £ where £ is 
the energy density, it is not possible to determine pa¬ 
rameters such as the internal energy U or entropy S 
from measurements on a single (interacting) 2D cloud. 
Moreover, other variables such as Tan’s contact param¬ 
eter C require derivatives with respect to the interac¬ 
tion parameter, which can only be obtained through 
combining measurements of clouds prepared with differ¬ 
ent interaction strengths hi[kpa 2 D) where kp = \f2Trn 
is the Fermi wavevector and a 2 D is the 2D scattering 
length. As described in the main text, the Helmholtz 
free energy F is given hy F = F/{NEp) = p — ^p/Ep. 
Using the Tan relations for the contact density in 2D 
C = dF/{d\na 2 D) [HHH] one can find the average en¬ 
ergy per particle U/{NEp) = p/2 — C/kp. With this 


it is possible to determine further thermodynamic pa¬ 
rameters provided the experiment is repeated at multiple 
values of a 2 D- As we are required to differentiate with 
respect to a 2 D at fixed T and we only have data at four 
unique values of a 2 D we necessarily introduce additional 
uncertainties in the derived quantities when compared to 
parameters such as chemical potential and temperature 
which can be obtained from a single cloud. This is es¬ 
pecially important for the extreme data sets and such 
measurements could be improved with additional data 
covering a larger range of a 2 D- However, as the overall 
contribution of the contact is small for the range of in¬ 
teraction strengths covered in these experiments we still 
find relatively small error bars and good agreement with 
the virial expansion at high temperatures. 

Having obtained U/{NEp) it is straightforward to 
extract the entropy per particle S/{NkB) from the 
Helmholtz free energy. In Fig. |S5[ b) we plot the entropy 
in this way. It is interesting that the measured entropy 
differs only very slightly from the ideal gas entropy for 
all interaction strengths we consider. Deviations from the 
ideal gas entropy may become more significant at lower 
temperature, particularly below the superfluid transition. 


[1] P. Dyke, K. Fenech, T. Peppier, M. G. Lingham, S. 
Hoinka, W. Zhang, B. Mulkerin, H. Hu, X.-J. Liu, and 
C. J. Vale, Phys. Rev. A 93, 011603(R) (2016). 

[2] G. Reinaudi, T. Lahaye, Z. Wang, and D. Guery-Odelin, 
Opt. Lett. 32, 3143 (2007). 

[3] L. Chomaz, L. Gorman, T. Yefsah, R. Desbuquois, and 
J. Dalibard, New J. Phys. 14, 055001 (2012). 

[4] T. Yefsah, R. Desbouqouis, L. Chomaz, K. J. Gunter, 
and J. Dalibard, Phys. Rev. Lett 107, 130401 (2011). 

[5] X.-J. Liu, H. Hu, and P. D. Drummond, Phys. Rev. B 
82, 054524 (2010/ 

[6] D. S. Petrov, and G. V. Shlyapnikov, Phys. Rev. A 64, 
012706 (2001). 

[71 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 










































5 



Figure S5. (a) Chemical potential and (b) entropy per particle of a two-dimensional Fermi gas with attractive interactions. The 
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for the corresponding interaction strengths. 
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